Topological Nature of the Phonon Hall Effect 
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We provide a topological understanding of the phonon Hall effect in dielectrics with Raman spin- 
phonon coupling. A general expression for phonon Hall conductivity is obtained in terms of the 
Berry curvature of band structures. We find a nonmonotonic behavior of phonon Hall conductivity 
as a function of the magnetic field. Moreover, we observe a phase transition in the phonon Hall 
effect, which corresponds to the sudden change of band topology, characterized by the altering of 
integer Chern numbers. This can be explained by touching and splitting of phonon bands. 
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Recent years have witnessed a rapid development of an 
emerging field - phononics, the science and technology 
of controlling heat flow and processing information with 
phonons Indeed, in parallel with electronics, vari- 
ous functional thermal devices such as thermal diode Q , 
thermal transistor Q , thermal logic gates 0] and thermal 
memory [f|, etc., have been proposed to manipulate and 
control phonons, the carrier of heat energy and informa- 
tion. However, different from electrons, phonons as neu- 
tral quasiparticles, cannot directly couple to the magnetic 
field through the Lorentz force. Therefore, it is a surprise 
that Strohm, Rikken, and Wyder observed the phonon 
Hall effect (PHE) - the appearance of a temperature dif- 
ference in the direction perpendicular to both the applied 
magnetic field and the heat current flowing through an 
ionic paramagnetic dielectric sample [6| . It was confirmed 
later by Inyushkin and Taldenkov [7J . Since then, several 
theoretical explanations have been proposed 0^3] to un- 
derstand this novel phenomenon. 

For electronic transport properties in various quan- 
tum, spin, or anomalous Hall effects ll-13j], topological 
Berry phase has been successfully used to understand the 
underlying mechanism 14] . Such an elegant connection 
between mathematics and physics provides a broad and 
deep understanding of basic material properties. How- 
ever, because of the very different nature of electrons 
and phonons, a topological picture related to the PHE 
is not straightforward and obvious, and therefore, is still 
lacking. 

In this Letter, we explore the topology of phonon 
bands in a two-dimensional honeycomb lattice with Ra- 
man type spin-phonon interaction. A general expression 
for phonon Hall conductivity in terms of Berry curva- 
ture is derived. The phonon Hall effect is not quantized, 
although the Chern numbers are quantized to integers. 
We find that there exists a phase transition associated 
with the PHE, due to the discontinuous jump of Chern 
numbers. 

We start with a Hamiltonian for an ionic crystal lattice 
in a uniform external magnetic field [l[ , which reads in a 



compact form as 
1 



H 



„(P _ Au) T (p - Au) + -u T Ku 
\p T P + \ uT ( K ~ ^ 2 ) M + u T Ap. 



(1) 



Here, u is a column vector of displacements from lattice 
equilibrium positions for all the degrees of freedom, mul- 
tiplied by the square root of mass, p is the conjugate mo- 
mentum vector, and K is the force constant matrix. The 
superscript T stands for the matrix transpose. A is an 
antisymmetric real matrix, which is block diagonal with 
h 



elements A = 



-h 



(in two dimensions), where h is 



proportional to the magnitude of the applied magnetic 
field, and has the dimension of frequency. For simplicity, 
we will call h the magnetic field later. The on-site term, 
u T Ap, can be interpreted as the Raman (or spin-phonon) 
interaction [16J. The Hamiltonian (1) is positive definite. 

By applying Bloch's theorem, we can describe the sys- 
tem by the polarization vector x — (/i, e) T , where /i and 
e are associated with the momenta and coordinates, re- 
spectively. The equation of motion can be expressed as 



d 
l dt X 



H e sX, H c 



-A -D 
I -A 



(2) 



where D(k) = -A 2 + J2i> iiz,;<e l(Ri ' -Ri) ' k is the dynamic 
matrix as a function of wave vector k; Kiji is the subma- 
trix between unit cell I and I' in the full spring constant 
matrix K; R/ is the real-space lattice vector; A is block 
diagonal with elements A, and I is an identity matrix. 
Here, D, A, Kij/, and I are all 4 x 4 matrices for the two- 
dimensional honeycomb lattice. The eigenvalue problem 
of the equation of motion @ reads: 



(3) 



where x a = (n ai e a ) T is the right eigenvector of the a- 
th branch and w a is the corresponding eigenfrequency. 
Because of the non-Hcrmitian nature of i/eff, the left 



2 



eigenvector is different, and is given by — (/t CT ,e CT ) = 
(ej., — /xj.) I [— 2iu) a ). The orthonormal condition is e^Cu' + 

By taking into account only positive eigenfrequency 
modes, displacement and momentum operators can be 
written in the second quantization form. From the def- 
inition of energy current density J = — 
Ri')uf Rimini [H H, H|, the current density vector can 
be expressed as 



J = Ji(a T a) + 3-2{a) a\aa) 



Here, Ji = w E 



k dk £fc ' a fe a fc' e 



: .(ui k —L 



(4) 



and J 2 = ^ £ 



fe,fe' 



w* ( £ fc dk ' e k' a 'k a 'k' e 



) -et f aiaie i ^"+ u k')t + 



el d -^e k ,a k a k ,e-^^)5^, where k = (k, a) 
considers both the wavevector and the phonon branch. 
It should be noted that the aW and aa terms also 
contribute to the off-diagonal elements of thermal 
conductivity tensor, although they have no contribution 
to the average heat flux. The diagonal term e\. dD J^ (k 

9k 



in 3 ± corresponds to w CT %f • Only the off-diagonal terms 



in Ji and J 2 contribute to the Hall conductivity, which 
can be regarded as the contribution from anomalous 
velocities similar to the one in the intrinsic anoma- 
lous Hall effect [IH • Using the Green-Kubo formula 

K *y = W Io h(iX ir dt ( JX (- iX ) JV V)) eq 0> one can 
obtain phonon Hall conductivity as 16]: 



8VT V 



f{u a )(u a +uv) 2 x 



e ^lfr eCT ' e a'l£' e<T (kx ^ k y ) 



, (5) 



where f(u) a ) — (e^/^B 7 )- l)" 1 is the Bose distribution 
function, V is the total volume of the sample, and the 
phonon branch index a here includes both the positive 
and negative eigenvalues without restrictions. It can be 
proved that the phonon Hall conductivity K xy satisfies 
the Onsager reciprocal relations (la) . 

In Fig. [T] we show the phonon Hall conductivity of hon- 
eycomb lattices calculated from Eq. ([3]). The parame- 
ters used in our numerical calculations are the same as 
in Ref. jij]. The coupling matrix between two sites is 
configured such that the longitudinal spring constant is 
Kl = 0.144eV/(uA 2 ) and the transverse one Kt is 4 
times smaller. The unit cell lattice vectors are (a, 0) and 
(a/2, a-s/3/2) with a = 1 A. 

It is found that when h is small, n xy is proportional to 
h [la) , while the dependence becomes nonlinear when h 
is large. As h is further increased, n xy increases before 
it reaches a maximum at certain value of h. Then n xy 
decreases and goes to zero at very large h. This can be 
understood as follows: numerical calculation shows that 
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FIG. 1: (color online) (a) Phonon Hall conductivity vs mag- 
netic field for different temperatures. The dotted, dashed, 
and solid lines correspond to T = 50, 100, and 300 K, respec- 
tively. The inset shows /i-dependence of K, xy at low tempera- 
tures: T = 10 (solid line), 20 (dashed line), and 40 K (dotted 
line), (b) d,K xy /dh as a function of h at different tempera- 
tures: T — 50 (dotted line), 100 (dashed line), and 300 K 
(solid line); here Nl = 400. The inset in (b) shows the h- 
dependence of dn xy /dh for different size Nl at T = 50 K, 
around h ~ 25.5 rad/ps; from top to bottom, Nl = 80,320, 
and 1280, respectively. 



u)„ ~ ah, which can also be obtained from the equation 
[(— iuj a + A) 2 + D~\e a = [la ], thus we can obtain approx- 
imately K xy ~ h 2 /(e l3nah - 1) from Eq. ©. In the weak 
magnetic field limit n xy cx h, while in the strong field 
limit, n xy — » 0. The on-site term A 2 in the Hamiltonian 
([l} increases with h quadratically so as to blockade the 
phonon transport, which competes with the spin-phonon 
interaction. Therefore, as h increases, n xy first increases, 
then decreases and tends to zero at last. At low temper- 
atures, K xy oscillates around zero with the variation of h, 
as shown in the inset in Fig. Ufa). 

There is a subtle singularity near h ~ 25 rad/ps in 
Fig.flja); we thus plot the first derivative of n xy with re- 
spect to h at different temperatures in Fig.[IJb). It shows 
that, at the relatively high temperatures, the first deriva- 
tive of phonon Hall conductivity has a minimum at the 
magnetic field h c ~ 25.4778 rad/ps for the finite-size sam- 
ple N L = 400 (the sample has N — N 2 unit cells). The 
first derivative dn xy /dh at the point h c diverges when the 
system size increases to infinity. The inset in Fig. (Ub) 
shows the finite-size effect. At the point h c , the second 
derivative d 2 n xy / dh 2 is discontinuous. Therefore, h c is 
a critical point for the PHE, across which a phase tran- 
sition occurs. At low temperatures, the divergence of 
dn xy /dh is not so evident as that at high temperatures. 
However, if the sample size becomes larger, the discon- 
tinuity of d 2 n xy /dh 2 is more obvious, as illustrated in 
Fig. [lib). For different temperatures, the phase transi- 
tion occurs at exactly the same critical value h c , which 
strongly suggests that the phase transition of the PHE is 
related to the topology of the phonon band structure. 



3 



In the following, we would like to connect the PHE 
with the Berry phase to examine the underlying topolog- 
ical mechanism. As is wellknown, the band structure of 
crystals provides a natural platform to investigate the ge- 
ometric phase effect. Since the wave-vector dependence 
of the polarization vectors is inherent to the Hall prob- 
lems, the Berry phase effects are intuitively expected for 
the PHE in the momentum space. Following Berry's ap- 

proach [14], we set x(t) = e l ^ {t) ' 1 lo dt '^ k ^ Xa (k(t)), 
and then insert it into Eq. ([2]). The Berry phase is ob- 



tained as 7 
Berry curvature emerges as 



§ A£ • dk, with A£ 



^ffllf ' an d the 



—AZ 



—AZ 



(6) 



where, 



n 



" dk x V' dk v t<T e ff dk y t<T fc (T' dk x to 



is the contribution to the Berry curvature of the band 
a from a different band a'. The associated topological 
Chern number is obtained through integrating the Berry 
curvature over the first Brillouin zone as 



2tt 



dk x dk y nz k =^Yjii k „, (8) 



BZ 



L 2 

k 



where, L is the length of the sample. The phonon Hall 
conductivity formula, Eq. ([5]), is recasted into 



8VT ^ r 



(9) 



Here V = L 2 a. The term (oj 



>) 2 relating to the 



phonon energy is an analog of the electrical charge term 
e 2 in the electron Hall effect, thus the phonon Hall con- 
ductivity Eq. © is similar to but different from the elec- 
tron case because the phonon energy term can not be 
moved out from the summation. Although the formula 
is derived from the phonon transport in the crystal-lattice 
system, we note that the thermal Hall conductivity for 
the magnon Hall effect [l9[ can also be cast into the form 
of Eq. ^ with a different expression for the Berry cur- 
vature. Therefore, the Hall conductivity formula can 
be universally applicable to the thermal Hall effect in 
phonon and magnon systems without restriction for spe- 
cial lattice structures. 

Without the Raman spin-phonon interaction, namely, 
h = 0, ^ rT k is zero everywhere and the phonon Hall 
conductivity vanishes. When a magnetic field is applied, 
the Berry curvature is nonzero, and consequently, the 
PHE appears. It is found that if the system exhibits 
symmetry satisfying SDS^ 1 = D, SAS^ 1 = —A (e.g., 
mirror reflection symmetry) , the phonon Hall conductiv- 
ity is zero 0, [HI ■ This symmetry principle can also be 
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CO FIG. 2: (color online) (a)-(d) The contour map of Berry 



curvatures for bands 1 — 4 at h c - = h c — 10~ rad/ps; (e)- 

(h) The contour map of Berry curvatures for bands 1 — 4 at 
h c + = h c + 10 -2 rad/ps. For (a)-(h), the horizontal and ver- 
tical axes correspond to wave vector k x and k y , respectively. 

(i) Q at different magnetic fields. The solid and dashed lines 
correspond to Q, 2 and il 3 at h c - respectively, while dotted 
and dash-dotted lines correspond to those at h c +. (j) Chern 
numbers of four bands: C 1 (solid line), C 2 (dashed line), C 3 
(dotted line), and C 4 (dash-dotted line), (k) The dispersion 
relation of band 2 and 3 at different magnetic fields in the 
vicinity of h c . The dashed, solid and dotted lines correspond 
to the bands at h c -, h c , and h c +, respectively. The lower 
three and upper three correspond to bands 2 and 3, respec- 
tively. k y — in (i) and (k). 



applied to the topological property of the phonon bands: 
we find that = provided that such symmetry ex- 

ists, such as in the square lattice system. Whereas if such 
symmetry is broken for the dynamic matrix, the system 
can possess nontrivial Berry curvatures. In the system 
with the PHE, if the magnetic field changes, the Berry 
curvatures are quite different. However, we find that the 
associated topological Chern numbers remain constant 
integers with occasional jumps when h is varied. There- 
fore, the Chern numbers given by Eq. ([5]) are topological 
invariant, which indeed illustrates the nontrivial topol- 
ogy of the phonon band structures. Although the Chern 
numbers are quantized to integers, the phonon Hall con- 
ductivity is not, due to the extra term /(w cr )(w (T + ui a ') 2 - 
Thus, the analogy to the quantum Hall effect is incom- 
plete. 

In the vicinity of the critical magnetic field h Cl we find 
that the phase transition is indeed related to the abrupt 
change of the topology of band structures. The Berry 
curvatures for different bands near the critical magnetic 
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field are illustrated in Fig. [2^a-h) . We find that with 
an infinitesimal change of magnetic field around h c , the 
Berry curvatures around the V (k = 0) point of bands 2 
and 3 are quite different, whereas those of band 1 and 4 
remain unchanged. To illustrate the change of the Berry 
curvatures clearly, we plot the cross section of the Berry 
curvatures along the k x direction for bands 2 and 3 in 
Fig. EJi), which shows explicitly that the Berry curva- 
tures change dramatically above and below the critical 
magnetic field h c . Below the critical point, the Berry 
curvature for band 2 in the vicinity of T point contributes 
Berry phase 2ir (—2ir for band 3), which cancels that from 
if, and K' points, so that the Chern number is zero for 
bands 2 and 3, as indicated in Fig.[2]j). However, above 
the critical point, the sum of Berry curvature at T point 
is zero, and only the monopole at if, and K' points con- 
tributes to Berry phase (— 27r for band 2 and 2tt for band 
3). Therefore, the Chern numbers jump from to ±1, as 
shown in Fig. EJj). This jump indicates that the topol- 
ogy of the two bands suddenly changes at the critical 
magnetic field, which is responsible for the phase transi- 
tion. From a calculation on the kagome lattice, which has 
been used to model many real materials [20l ] , we also find 
qualitatively similar phase transitions due to the sudden 
change of topology, where the phonon Hall conductivity 
has three singularities of divergent first derivatives cor- 
responding to three jumps of the Chern numbers. 

To further investigate the mechanism of the abrupt 
change of the phonon band topology, we study the dis- 
persion relation near the critical magnetic field. From 
Fig. HJk), we can see that band 2 and 3 are going to 
touch with each other at the T point if the magnetic 
field increases to h c ; at the critical magnetic field, the 
degeneracy occurs and the two bands possess the cone 
shape; above the critical point h c , the two bands split 
up. Therefore, the difference between the two bands de- 
creases below and increases above the critical point h c . 
The property of the dispersion relation in the vicinity 
of the critical magnetic field directly affects the Berry 
curvature of the corresponding bands. 

In summary, we have studied the PHE from a topolog- 
ical point of view. By looking at the phases of the polar- 
ization vectors of both the displacements and conjugate 
momenta as a function of the wave vector, a Berry curva- 
ture can be defined uniquely for each band. This Berry 
curvature can be used to calculate the phonon Hall con- 
ductivity. Because of the nature of phonons, the phonon 
Hall conductivity, which is not directly proportional to 
the Chern number, is not quantized. However, the quan- 
tization effect, in the sense of discontinuous jumps in 
Chern numbers, manifests itself in the phonon Hall con- 
ductivity as a singularity of the first derivative with re- 
spect to the magnetic field. 

The topological approach for phonon Hall conductivity 
proposed here is general and can be applied to the real 
materials in low temperatures where the thermal trans- 



port is ballistic. It can also b e ap plied to the magnon 
Hall effect discovered recently [19(. Phase transition in 
the PHE, explained from topological nature and disper- 
sion relations, can also be generalized to study the phase 
transition in other Hall effects and/or noncquilibrium 
transport. In line with recently reported Berry-phase- 
induced heat pumping [2~H and the Berry-phase contri- 
bution of molecular vibrational instability [22|, we hope 
our present results do invigorate the studies aimed at 
uncovering intriguing Berry phase effects and topological 
properties in phonon transport, which will enrich further 
the discipline of phononics. 
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SUPPLEMENTARY INFORMATION FOR "TOPOLOGICAL NATURE OF PHONON HALL EFFECT 



In this supplementary material, we discuss the origination of the Hamiltonian [Eq. (1) in the main text] in the first 
section; then we present the detailed derivation of the general formula of phonon Hall conductivity in terms of Berry 
curvature, in which we also give the explicit expression for the dynamic matrix and give the proof for the symmetry 
principle. Finally, we discuss the numerical calculation for Chern numbers. 



In the presence of a magnetic field, according to Ref. the kinetic energy of each site in ionic crystal lattices 
without free charges is expressed as: 



where, r a — R Q + u Q /^/m a , R Q is the equilibrium coordinate of the ion at site a, and u Q denotes the displacement 
multiplied by the square root of the ion mass m a . p a is the corresponding momentum divided by the square root of 
mass m a . q a is the ionic charge at site a. A a denotes the electromagnetic vector potential, which, using the Lorenz 
gauge condition, can be related to the ionic displacement as [l[ 



ABSTRACT 



DISCUSSION ON THE HAMILTONIAN 




(SI) 




(S2) 



Thus, Eq. ([Sljl is recasted as: 



T a = 




2m, 



a. 



B x u a | 2 . 



(S3) 



If the magnetic field with magnitude B is applied along z direction and we only consider the two-dimensional (x and y 
direction) motion of the system, then the kinetic energy of ion a can be expressed (it is straightforward to generalize 
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to high dimentions) as: 

T a = i(p a - ^ a u a ) T {p a - A a u a ), (S4) 

where p a = (p a x,Pay) T , u a = (u ax ,u ay ) T , and A a = ^ ^ '(T ^ ' wnere ^ a ~ ~~ QaB /(2m a ). Note that there are 

both positive and negative ions in one unit cell. For a general ionic paramagnetic dielectric, mostly, the mass of the 
positive ion is larger than that of the negative one. For instance, in the experimental sample TbaGasO]^, the ratio 
m(+q)/m(—q) is about 4.3 in one unit cell. Therefore the negative ions will dominate in the contribution to h a , which 
makes h a have the same sign as that of the applied magnetic field B. Under the mean-field approximation, we can 
set h a = h, which is site-independent and is proportional to the magnitude of the applied magnetic field. 

Combining the kinetic energy with the harmonic inter-potential energy, we can write the whole Hamiltonian as 

H=^{p- Au) T {p - Au) + ^u T Ku 7 (S5) 

where A is an antisymmetric real matrix with block-diagonal elements A a . u and p are column vectors denoting 
displacements and momenta respectively, for all the degrees of freedom. K indicates the force constant matrix. 
Finally, after the rearrangement, we have 

H = ^p T p+^u T {K - A 2 )u + u T Ap, (S6) 

which is exactly the second row of Eq. (1) in the text. 

The Hamiltonian Eq. (|S6[) [Eq.(l) in the main text] is essentially the same as that used in Ref. resulting from 
the phenomenological Raman interaction. The only difference is the term proportional to A 2 which makes the above 
Hamiltonian positive definite. The Raman interaction, proposed to study spin-phonon interactions (SPI) based on 
quantum theory and fundamental symmetries 0^], can be expressed as 

ff/ = js-(uxp). (S7) 

Here, g denotes a positive coupling constant, and s is the isospin for the lowest quasidoublet. In the presence of 
a magnetic field B, each site has a magnetization M. For isotropic SPI, the isospin s is parallel to M, and the 
ensemble average of the isospin is proportional to the magnetization, which can be expressed as (s) = cM with c the 
proportionality coefficient (Ref. [2-5]). In the mean- field approximation, the Raman type SPI reduces to 

ff/ = h-(uxp), (S8) 

where h = gcM, and M is proportional to the magnetic field B. If the magnetic field is applied along the z direction, 
then the SPI can be written as 

H I = u T Ap. (S9) 

By treating the phonon system under harmonic approximation, the total Hamiltonian for the whole lattice can be 
written as (Ref. 

H = \p T P+ ^u T Ku + u T Ap. (S10) 

Note that this Hamiltonian Eq. (|S 10[) is not positive definite. In Ref. [J], the authors added an arbitrary onsitc 
potential in order to make the Hamiltonian positive definite. However, in the calculation of phonon Hall effect for 
the four-terminal junctions, such non-positive-definite Hamiltonian does not cause any problem because the thermal 
junctions will stabilize the system 

From the first physical picture of spin-phonon interaction in ionic crystal lattice with an applied magnetic field 
(Eq. ??^ IS6I) . the additional term proportional to A 2 emerges naturally to make the Hamiltonian positive definite. 
Therefore, in this work we choose the positive definite Hamiltonian Eq. (|S6p [Eq. (1) in the main text]. 
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PHONON HALL CONDUCTIVITY FROM GREEN-KUBO FORMULA 

The Hamiltonian Eq. (|S6j) is quadratic in u and p, and we can write the equation of motion as 

p = — {K — A 2 )u — Ap, (Sll) 
u = p-Au. (S12) 

The equation of motion for the coordinate is, 

u+ 2Au + A 2 u + (K - A 2 )u = 0. (S13) 
Since the lattice is periodic, we can apply the Bloch's theorem it/ = ee i ( R i k - wt ). The polarization vector e satisfies 

[(-iuj + A) 2 + D]e = 0, (S14) 

where D(k) = —A 2 + J^>, Ku' e l ( R i' k denotes the dynamic matrix and A is block diagonal with elements A. 
D, Ki p, and A are all nd x rod matrices, where ro is the number of particles in one unit cell and d is the dimension of 
the motion. 

To calculate the dynamic matrix -D(k), we give an example for the two-dimensional honeycomb lattice, where 
11 = 2, d = 2. We only consider the nearest neighbor interaction. The spring constant matrix along x direction is 

<S15) 

Kl = 0.144eV/(uA 2 ) is the longitudinal spring constant and the transverse one Kt is 4 times smaller. The unit cell 
lattice vectors are (a, 0) and (a/2, ay/Z/2) with a = 1 A. 

To obtain the explicit formula for the dynamic matrix, we first define a rotation operator in two dimensions as: 

Uos 6 -sin fly 
' V sine' cos# J 

The three kinds of spring-constant matrices between two atoms are K$i — U(tt/2)K x II(—tt/2), K02 = 
U(tt/Q)K x U(—tt/6), K03 — U(~tt/6)K x U(tt/6), which are 2x2 matrices. Then we can obtain the on-site spring- 
constant matrix and the four spring-constant matrices between the unit cell and its four nearest neighbors as: 

jy- ( ^01+^02 + ^03 —if 02 \ fQta\ 

K °-{ -K 02 K al +K a2 +K a3 )> ^ 16 > 



*-U.S)- K >-{-lD- "-(l-?)- Moi'" 1 )' (S17 » 

which are 4x4 matrices. Finally we can obtain the 4x4 dynamic matrix D(k) as 

£>(k) = -A 2 + K + K ie lk * + K2e ^/2+V3k y /2) + Kae -ik* + X4e - J (fc,/2+v^fe y /2) j (g 18 ) 

where, A 2 = —h 2 ■ I, and I denotes the 4x4 identity matrix. 

Equation (|S14I) can be written in a form as a standard eigen-problem given in Eq. (3) in the main text, if we rewrite 
the equations of motion. Using Bloch theorem, Eqs. (IS11[) and (|S 12[) can be recasted as: 

i-^x = H eS x, H eS = i ( f ^ • ( S19 ) 

x = (/i, e) T is the polarization vector , where column vectors p, and e are associated with the momenta and coordinates 
respectively. Therefore, the right eigenvector and left eigenvector satisfy: 

H eS x a = u a x<j, x^Heff = io a x T a . (S20) 
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where the right eigenvector x a = (^ cr ,e cr ) T , the left eigenvector — ( e t 5 _ )/(— 2iu; CT ), and a indicates the branch 
index. Because the effective Hamiltonian H e g is not hermitian, the orthonormal condition then holds between the 
left and right eigenvectors. The eigenmodes can be normalized as x^x a = 1, which is equivalent to [IJ 



4 ta + —elAe a = 1. (S21) 
To solve the eigensystem, we require the following relations: 

e -k,-er = e k,o-5 ^-k,-cr = — Wk,cr. (S22) 

In the following, we use k = (k, a) to specify both the wavevector and the phonon branch. By taking into account 
only positive eigen-modes (w > 0), displacement and momentum operators are taken in the second quantization form: 



u t = E^ el(Rrk ^ fct) v^^+ h - c -' ( S23 ) 



where a > 0, a^ is the annihilation operator, and h.c. stands for hermitian conjugate. The momentum and displace- 
ment polarization vectors are related through fi^ = —iuk^k + ^.efc- We can verify that the canonical commutation 
relations are satisfied: [ui,pf,] = ifiSij/I, and H = TiuJk(a\.ak + 1/2). 
The energy current density is defined as [9(: 

3 = ^E( R '- R '')"F^M'«Z'> ( S25 ) 

1,1' 

where V is the total volume of N unit cells. The current density vector can be expressed in terms of the cre- 
ation/annihilation operators as 

J = Ji (at a) + J2(ftV, aa); 

Ji = w (v^ + v^) 4^*4^^^^—^; (S26) 



fe.fe' 

We note that the a^a^ and aa terms also contribute to the off-diagonal elements of the thermal conductivity tensor, 
although they have no contribution to the average energy current. Based on the expression of heat current, the 
phonon Hall conductivity can be obtained through the Green- Kubo formula [loj : 



h/(k B T) 

d\ I dt(J*{-i\)JV(t)\, (S27) 



_ V_ 

Kxy ~hT 

where the average is taken over the equilibrium ensemble with Hamiltonian H. Substituting the expression J into 
Eq. (|S27p . the phonon Hall conductivity is obtained as 

K , „(*). 

= m Io /(kBT) f °°dt < Jf (-iX) J?(t)> eq ; (S28) 



*® = wrIo /ikBT) d\J °°dt(j%(-i\)JUt)) 



eq 



Note that the averages of the cross terms (jf (— iA) J* (t)} and ( Jf (— iX)J\ (t)) e are zero. 
First we calculate the term . Combining the result 

( a l a j a l a i)cq = fifkSijSki + fi(fj + l)8uSjk, (S29) 
where /j = [eP" 1 " 1 — 1) _1 is the Bose distribution function, with the result 
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which is obtained by differentiating the Eq. (S20) (Eq. (3) in the main text) and x^x a = 1, we obtain 

(1) n ST Iff \ fl Ml/ i M2 1 a gfcg g CT ' 8fc y J zoom 

k,it>0,<t'>0 

Because of Eq. (|S22I) and the following property: 

D ab (-k) = D* 6 (k) - D ba (k), (S31) 
we can transform from the positive-frequency bands to the negative-frequency band, and obtain 

* ■ f t JLD , f t ao. 



We can prove = —K yx , such that 



Kxy 16VT 



ft M f ,ft l£ f _ft M f ,ft 3D 



Because of = — £1% X i the phonon Hall conductivity can be written eventually as 

h 



(S33) 



Here, a, a' can be both positive or negative. 
Second, we calculate n ab . Utilizing the results 

{a\a)a k ai) e ^ = fifj(S lk Sji +S u 8 jk ); 

( a i a j a Wl)cq = (1 + + fj)(8ik8jl + SuSjk), 

and the relation /(—a;) = —1 — /(w), after some algebraic derivation similar to the above, we obtain 

+ ■ e t9fl £ , e t ao 

k.crcr'<0 

Therefore, the total phonon Hall conductivity can be written as 

n [W \ \i / . n2 * ° 9fex t,T V dk y , S 

= Wf2Z - + 4w K-av) 2 • (S35) 

V rr=trr' v 7 



^ [/K)-/K0]K + ^) 2 ^X> (S36) 



fox f$V A , . , . ( , > , ■ . \2 



^ /k)k+^) 2 ^:1- b . (S38) 



where is the total volume of N = N\ unit cells. In the above formula, the phonon branch a includes both positive 
and negative values without restriction. We start with the positive frequency bands to derive the conductivity formula. 
Through some transformations, we finally obtain the simplified formula for phonon Hall conductivity which combines 
the contribution from all the frequency bands. The formula Eq. (|S38[) is different from that given in Ref. In 
Ref. Q the contribution for phonon Hall conductivity from J 2 was omitted, which is incorrect. 
From the Eq. (|S14I) , we obtain 

e -k,o-(-^) = £k,<r(A); w_ k , CT = Wk.cr, (S39) 
and because of D(k) = D*(-k), e£§j£e£, = e t>mt e °-> we have 

fll a ' k %-A) = il a k a ' k (~k,A) = -Q^' k (— k, A). (S40) 
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FIG. SI: (color online) (a) Phonon Hall conductivity vs applied magnetic field for two-dimensional honeycomb lattice, (b) 
Phonon Hall conductivity vs temperature at fixed magnetic field h = 1 rad/ps. The inset of (b) shows the product of phonon 
Hall conductivity and temperature K xy T vs magnetic field h for different temperatures. 



So we obtain 

Kxyi-A) = K yx {A) = —K xy (A). (S41) 

The Onsager reciprocal relations are satisfied. 

If the system possesses the symmetry which satisfies 

SDS- 1 = D, SAS- 1 = -A, (S42) 

where S represents any symmetric operation, and from Eq. (|S14[) . we obtain 

Se(A) = e(-A). (S43) 

Using the definition of the dynamic matrix D = —A 2 + iY/ j ;<e J ( Ri '~ R '* l ' k and SDS- 1 = D, we can obtain 

S W a S - 1 = W Q > (* = X ^ (844) 
Inserting S~ 1 S = I into Eq. (|S37|) . we obtain 

^Z'k y {~ A ) = ^Z'k v (A) ■ (S45) 

Then it is easy to obtain K xy (—A) = K xy (A), and because of the Onsager relation, one can easily obtain that 

K xy = 0, if SDS' 1 = D, SAS' 1 = -A. (S46) 

Fig.lSlTa) shows the phonon Hall conductivity with magnetic field for different temperatures. In the weak magnetic 
field range, the phonon Hall conductivity k xy is proportional to the magnetic field, which is consistent with all 
the experimental and theoretical results. We plot phonon Hall conductivity with a large range of temperatures in 
Fig. ISir b). At very low temperatures, the phonon Hall conductivity is proportional to 1/T. k xy T will be constant 

(2) 

for different temperatures lower than IK. This is due to the contribution from n xy : if T — > 0, 1 + / — > 1, then 
the conductivity linear with 1/T tends infinity. While the longitudinal thermal conductivity k xx is infinite for any 
temperature thus when T — > 0, the transverse Hall conductivity, n xy — > oo, has the ballistic property similar to the 
longitudinal one. If temperature is very high, all the modes contribute to the thermal transport, and / ~ ksT/iJvj)), 
then the phonon Hall conductivity becomes a constant, which can be seen in Fig. ISllf b) . 
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THE BERRY PHASE AND BERRY CURVATURE 

Using the similar method proposed by Berry we derive the Berry phase and Berry curvature in the following. 
Starting from 

i^x{t) = H eS x(t) (S47) 

and substituting 

we can obtain the Berry phase across the Brillouin zone as 

la = j> A CT (k)dk, A ff (k) = iZa-^x*. 
Here x a ,x^ correspond to the right and left eigenvectors, and x^x a ' = , Ylx^x^ = I. A°"(k) is the so-called 



Berry vector potential. Therefore the Berry curvature is obtained through the Stokes theorem as: 

~ r T gHtft _ =T dH, 
L a dk x iff dk 



W k k = -^-Al - -£-A% = i Y — , dky ,„ (S48) 

k * k " dk x dk v k * (u a - uj a ,) 2 v ' 



Inserting the vector x and the expression of matrix H e g, we obtain 

El (7 r)k a ^o' rlk a ""HI X — T ' 
— — — - 9 = y nr k (S49) 

where Sl k a k indicates the contribution to the Berry curvature of the band a from a different band a' . Therefore, the 
phonon Hall conductivity formula Eq. (|S38[) can be interpreted in terms of the Berry curvature. 

THE CALCULATION OF THE CHERN NUMBER 

The topological Chern number is obtained by integrating the Berry curvature over the first Brillouin zone as 



C a = ^- { dk x dk y n a k k . (S50) 
/bz " " 



2tt 



For numerical calculation, we use 



L 2 

k 



where j^^l = j czttV an< ^ ^ = L 2 a, L 2 is the area of the sample. 

k 

To calculate the integer Chern numbers, large k-sampling points N is needed. However there is always a zero 
eigenvalue at the T point of the dispersion relation, which corresponds to a singularity of the Berry curvature. 
Therefore, we cannot sum up the Berry curvature very near this point to obtain Chern number of this band, unless 
we add a negligible on-site potential \u T V ons it c u to the original Hamiltonian. In Fig. IS2T a) . without the on-site 
potential, the Chern number of the fourth band is not an integer, no matter how large the sample size N = is (see 
Fig. IS2fb)). If we add the external on-site potential, the Chern number of the fourth band will become integer. In 
Fig. IS2f a). the C4 changes gradually to —1 with increasing the on-site potential, while other Chern numbers do not 
change. And from Fig. IS2f b). we see that with larger on-site potential, the Chern number of the fourth band could 
be an integer for smaller sample sizes. 

t Electronic address: renjie@nus.edu. sg 
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FIG. S2: (color online) (a) The four Chern number vs onsite potential V ona itc- The unit for onsite potential is longitudinal 
spring constant Kl- Here N = iV£ = 160000; (b) The Chern number of the fourth band changes with Nl for different onsite 
potentials. For both (a) and (b), h — 1 rad/ps. 
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